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Abstract. We solve the double exchange model on nanoscopic clusters exactly, and specifi- 
cally consider a six-site benzene-like nanocluster. This simple model is an ideal testbed for 
studying magnetism in nanoclusters and for validating approximations such as the dynam- 
ical mean field theory (DMFT). Non-local correlations arise between neighboring localized 
spins due to the Hund's rule coupling, favoring a short-range magnetic order of ferro- or 
antiferromagnetic type. For a geometry with more neighboring sites or a sufficiently strong 
hybridization between leads and the nanocluster, these non-local correlations are less rele- 
vant, and DMFT can be applied reliably. 
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1 Introduction 

In recent years, manganites such as Lai-^Ca^MiiOs 
(LCMO) have attracted great interest, mainly due 
to the colossal magnctorcsistancc [T] . In these ma- 
terials, the crystal field splits the five 3d orbitals 
of the Manganese atoms into two e g and three t2 g 
orbitals due to the perovskite structure. The latter 
are localized and as a consequence of Hund's ex- 
change half-filled, forming a spin 3/2. This spin, the 
itinerant e g electrons and their coupling, again by 
Hund's exchange, constitute the double exchange 
or ferromagnetic Kondo lattice model [2J. This ar- 
guably simplest model for manganites gives rise, in 
the bulk, to a ferromagnetic double exchange since 
a hopping energy of the e g electrons can be gained 
only for a ferromagnetic (FM) alignment of the tig 
spins. For half-filled e g bands one the other hand, 
the alignment of the ti g spins is antiferromagnetic 
(AF) due to superexchange. 

More recently, nanoclusters of manganites have 
been synthesized; and remarkably, a size-control 
of the charge and magnetic ordering in half-doped 
LCMO has been demonstrated [cOS]. Hence, size 
can be utilized to optimize the magnetic proper- 
ties and the magnetoresistance of manganites for 
technological applications. 

Not only these experiments, but also generally 
the emergence and peculiarity of magnetism in nan- 



oclusters, motivates us to study the double exchange 
model: 



H = —t 



2j]T Si S i . 



G) 



Here, c\ G (c^) are the creation (annihilation) op- 
erators of an electron at site i with spin a, and t 
is the effective hopping amplitude corresponding to 
the double exchange process. S and s are the spins 
of the localized and mobile electrons respectively, J 
is the Hund's coupling between those spins. Hence- 
forth we restrict ourselves to a Ising symmetry of 
the localized spins so that Si simplifies to Si = ±1, 



and 



2 X/cr= 



±i 



Note the results only 



depend on \S\J so that taking Si = ±1 instead of 
Si = ±3/2 only corresponds to a redefinition of J. 

Besides these usual terms of the double exchange 
model, we consider a hybridization Vi^k connect- 
ing each site i of the nanostructure to some non- 
interacting environment rj with eigenenergies e n k 
where k labels the different states in the respective 
environment, see Fig.[T] Such term allows one to de- 
scribe, e.g., a surface or electrodes (leads) applied 
to the system for electric transport measurements. 
Altogether, this leads to the Hamiltonian 
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2 Method 




Fig. 1: Schematic representation of the model sys- 
tem consisting of six sites, each of them is con- 
nected to a lead via the hybridization strength 
Virik = V8i n , and electrons can hop between the 
sites i and j of the nanostructurc via hopping chan- 
nels tij. Itinerant electrons interact with localized 
spins S at each site via Hund's coupling J. 



In analogy to Refs. [3] we consider a nanoscopic 
system made of six-sites arranged in a ring struc- 
ture. For the sake of simplicity, we will restrict 
ourselves to a configuration in which all sites are 
equivalent and each site is connected to its own 
paramagnetic metallic lead, i.e., V^k = V5i v , and 
take a constant density of states p = 1/(2D) for 
the non-interacting bath, where D is the bath half- 
bandwidth. As we will see, the competition be- 
tween kinetic energy, hybridization, and Hund's ex- 
change exhibits some very interesting properties, 
and can lead, depending on the parameters, to fer- 
romagnetic (FM) or antiferromagnetic (AF) (and 
of course paramagnetic) short-range ordering of the 
localized spins. 

In this paper, we solve the model ([2]) exactly on 
a six-site ring, and take it as a benchmark for test- 
ing the validity of a dynamical mean field theory 
(DMFT) [B] approximation. Let us note that the 
DMFT solution for the Kondo lattice model in the 
bulk was derived before by Furukawa [7j. In order 
to apply DMFT to a finite system, we employ its 
recently introduced nanoscopic version [5], which 
has been shown to be suitable to deal with complex 
nanostructures. The outline of the paper is as fol- 
lows: In Sec. [5] we present the path integral method 
for the exact solution, and the DMFT scheme. In 
Sec. l2.ll we present the results obtained solving the 
model exactly, and we use them as a benchmark 
for DMFT in Sec. 12.21 Finally, we summarize our 
results in Sec. 



Let us start with some simple symmetry considera- 
tions. In a six-sites system, with two possible align- 
ments for each localized spin Si, there are 2 6 = 64 
different possible configurations {Si}. Since in the 
(high temperature) paramagnetic phase the Z2 ro- 
tational symmetry is conserved, the overall align- 
ment of the spins does not matter, thereby reduc- 
ing the number of inequivalent configurations by a 
factor of two. In the following, and 1 will repre- 
sent the two different orientations of the localized 
spins. So, e.g., a configuration with only one spin- 
up would be represented by {100000}, which, in the 
paramagnetic phase, is also equivalent to {011111}. 

In order to check the role of the connectivity, 
in the following we will consider similar cases as in 
Ref. [5]: (i) hopping is only possible to next neigh- 
bors (" NN t" ) in which the number of inequivalent 
spin configuration is further reduced to eight, by 
symmetry, and (ii) hopping is equally possible to 
all sites ("all t"), with only four inequivalent spin 
configurations. The latter hopping topology is not a 
realistic one, as the hopping amplitude is a decreas- 
ing function of the inter-site distance. However it 
will be useful to study, in the philosophy of Ref. [5] , 
the effect of enhanced connectivity without chang- 
ing the geometry of the nanostructure. All these 
symmetry considerations are summarized in Table 

m 

For an isolated nanostructure, i.e., without hy- 
bridization to the leads, the double exchange model 
can be solved exactly by calculating all eigenvalues 
and eigenstates of the Hamiltonian ([2]). This proce- 
dure becomes however computationally unfeasible 
if the number of sites of the nanostructure is very 
large, or if non-interacting leads are included, as 
it would require many additional bath-sites for an 
accurate description of the leads. We therefore em- 
ploy a Green's function technique, were the non- 
interacting degrees of freedom can be integrated 
out exactly. In the Grassmann path-integral for- 
mulation [5] the expectation value of an observable 
O is calculated as 

(0>4E [Oe^, (3) 

where Z = Ylis } f e ~ S ^ * denotes the partition 
function of the system and the integration sym- 

boi / = n iCT / T>\4a , c if ] ru, 1 v K fcf! % J ^- 

resents the functional integral extending over the 
Grassmann variables cj^ (c icr ) and ad fe(7 (a r)fccr ) as- 
sociated to the fermionic creation (annihilation) op- 
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NN t 


all t 


configuration 


degeneracy 


tt of AF bonds 


degeneracy 


tt of AF bonds 


{000000} 


2 


0/6 


2 


0/15 


{100000} 


12 


2/6 


12 


5/15 


{110000} 


12 


2/6 


30 


8/15 


{111000} 


6 


2/6 


20 


9/15 


{101000} 


12 


4/6 




{100100} 


6 


4/6 


{101100} 


12 


4/6 


{101010} 


2 


6/6 



Table 1: Non-equivalent localized spin configurations for the six-site ring for both NN t and all t hopping 
topologies, taking into account Zi rotational symmetry and translational hivariance. 



erators. Finally, the action is given by 
5 {S ' } =/^[^cL(r)(9 r - M )c l(T (r) 

+ E 4 a {r)t ijCj(7 {r)-2Y J Jsi{T)S i 

+ E a Ua ( T ) ( d r + Cjfc ~ ( T ) 

r/ka 

+ J2 V ^(4a(^%kA-r) + h.c.)], (4) 

irjha 

where r S [0,/3) is the imaginary time, /3 = 1/T is 
the inverse temperature, and /z is the equilibrium 
chemical potential. Hence the Green's function is 
defined as 



G mna (j) = 7? E / C ma( T ) c L(0) 



i S{Si} . (5) 



Since the non-interacting electron operators only 
enter quadratically in the action (0}, we can inte- 
grate them out by a simple Gaussian integral [5], 
yielding the effective action 

S { £ } = f dr\Y,ci{T){dr ~ H)c ia {r) 

+ E c l(r)UjC ja (r) - 2jJ2S iSi (T)) 

ija i 
f/3 r/9 



/ drf dr' E c L( t )Aj(t - r/)c v (r') 



Here the hybridization function Zi^ (t) takes into 
account the virtual processes of itinerant electrons 
hopping back and forth between the leads and the 
nanostructure. Its Fourier transform, analytically 
continued to the real axis, corresponds to the re- 
tarded function 



a» = e 



v ivk v*, k 



ijk 



u + 15 — e 



ijk 



(7) 



At this point it is convenient to recast the Green's 
function of the nanostructure in terms of the 
effective action ([6]) considering the contribution of 
each configuration {Si} of the localized spins ex- 
plicitly 

G m „.(r) = E^ {Sl} G^(r), (8) 

{Si} 

via a functional integral for each configuration {Si} 
G££(t) = / c ma (r)cU0)e- S ^' } (9) 

where we now only integrate the interacting degrees 
of freedom, i.e., / = \[ ia J V[c[ a ,c ia ). The weights 
of the configuration-dependent Green's functions in 
Eq. ([5]) are defined by 



P {s,} = 



E^ }: 

{Si} 



(10) 



where 



z = e^ } = e/ 

{S,} {S t } J 



,-s 



{Si} 



(11) 



The calculation of the probabilities Eq. (fit))) re- 

quires the evaluation of the determinant of S^ s . 
This is more conveniently done in Fourier space 
where it takes the matrix form, in site indices i 
and j, 



(6) S ij(r (u)) = (-w-n + Aij^w)- JaSi)Sij+tij. (12) 



Calculating these determinants (probabilities) as 
well as the Green's function for the inequivalent 
of the 2 6 spin configurations {Si} yields the exact 
summation of the Green's function of the nanoclus- 
ter coupled to non-interacting leads. 

Besides the exact solution, we have also em- 
ployed DMFT. There are several versions of DMFT 
for treating nanoscopic or spatially inhomogeneous 
systems [5U9 UTUUTT] with minor differences. Here, 
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When only the kinetic term of Hamiltonian ^ 
is taken into account, the nanostructurc is trans- 
lational invariant and the wavevectors k = rmr/3 
(m = 0, . . . , 5) are conserved quantum numbers. 
Setting the chemical potential p = 0, yields the dis- 
persion relation 



we perform our calculations in the same fashion as 
in Rcf. [5] , restricting ourselves to the 1-particle, i.e. 
DMFT, realization of the dynamical vertex approx- 
imation [12013] applied to nanoscopic systems. The 
starting point of the DMFT scheme is the Green's 
function of the nanostructurc, of which the generic 
matrix element of its inverse reads 

{ G_1 }yv( a; ) = w ^ ~ *« ~ A a^) - ^VM, 

(13) 

where £ denotes the self-energy. For each of the 
(incquivalcnt) sites of the nanostructure, we define 
a local problem determining the associated Weiss 
field Gf 1 = [Gu] 1 + Sa out of the local block of 
site i of the Green's function Gu, with an initial 
guess of the (local) self-energy Uu (usually zero). 
This procedure maps each atom onto an effective 
one-body local problem where the local spin can ei- 
ther point up or down. Let us note that the DMFT 
approximation includes the local spins. That is, 
DMFT substitutes the actual localized spin con- 
figuration on the other sites by an effective Weiss 
field Q° , which in the paramagnetic phase is spin- 
independent. Hence localized spin up and down have 
the same probability and the interacting DMFT 
Green function of site i reads |7| 

1/2 1/2 

Gi(w) = p- + P . (14) 

Qf >) + J fi? (w) - J 

Since this has the form of a non-interacting Green's 
function averaged over a potential ± J, Eq. (|T4)l is 
temperature independent (which does not hold for 
DMFT Green's function in general). Mathemati- 
cally this is equivalent to binary disorder. Hence, 
Eq. (TI"4"]) allows us to determine a local self-energy 

Su = Qf 1 + Gi^ 1 for each site i. Using the calcu- 
lated self-energies as an input for the Dyson equa- 
tion (|13p we can finally calculate a new Green's 
function for the nanostructurc, and iterate the cycle 
sclf-consistcntly until convergence. [5] In the case 
we are considering all sites are equivalent due to 
the symmetry of the problem. Therefore we need 
to solve only one local problem, yielding a local 
self-energy which is the same for each site of the 
nanostructure. 

2.1 Results: Exact Solution 

In this section the exact results obtained are pre- 
sented. As already mentioned, there are several com- 
peting energy scales in the problem. Therefore it is 
useful to describe the non- interacting (J/t = 0) 
and isolated (V/t = 0) "benzene" molecule first, 
and analyze separately the effects of the hybridiza- 
tion and of the interaction. The following analysis 
is summarized in the spectral function, connected 
to the retarded Green's function via the relation 
A{uj) = -~Gu(lj + i5), shown in Fig.[2J 



E(k) = -2b cos(fc) - 2b' cos(2fc) - b" cos(3fc), (15) 

where, b, b', and b" are the nearest neighbor, the 
next-nearest neighbor, and the next-next-nearest 
neighbor hopping amplitudes, respectively. In the 
following t sets our unit of energy, and the hop- 
ping configurations as defined by: (i) NN t case: 
H = b" = 0, and (ii) all t case: b' = b" = b. The 
spectral function associated to this system is con- 
stituted by (5-likc peaks corresponding to the six 
energy eigenstates (some of which are degenerate) , 
and it is shown in the left panel of Fig. [2] for the 
NN t case (solid line). It is important to note that 
in both the NN t and all t case the spectral weight 
at the Fermi energy (Ep) is vanishing and the sys- 
tem is a band insulator. 

In order to show what happens when the struc- 
ture is connected to metallic leads, we restrict our- 
selves to the NN t case, but analogous consider- 
ations apply to any hopping topology. As already 
mentioned, the leads are bulk systems of non- interacting 
electrons described by a flat density of states p = 
1/(213). We consider the broad-band limit D 3> b 
(in order to avoid the particular physics arising at 
the edge of the leads' conduction band) and we set 
D = 20t in the following. The most important ef- 
fect of the hybridization between the ring and the 
non-interacting environment is the broadening of 
all peaks in the spectral function which arc within 
the bandwidth of the leads, while the shift of the 
eigenenergies, i.e. the real part of the hybridization 
function, is 0(w/D). With increasing V, the dis- 
crete spectrum of the molecule evolves into a broad 
band (see Fig. [3J left panel). 

The next step is to study what is the effect of 
the Hund's coupling J on such a nanoscopic system. 
In general, for a system of N sites, there are 2^ dif- 
ferent configurations of the localized spins, which 
means that there are 2N2 N eigenstates. Most of 
them are degenerate in the non-interacting, para- 
magnetic system with translational invariance. In 
the presence of J and for a given configuration 
{Si}, translational invariance is broken, and some 
of those degeneracies are lifted. At low enough tem- 
peratures, however, only the spin configurations cor- 
responding to the lowest energy will be populated, 
so that the main effect of J is to increase the size of 
the gap. For large J, one can show that the size of 
the gap is 2 J (see Fig. [2j middle panel). At higher 
temperature (/?£ = 5) , due to the broadening of the 
Fermi-Dirac distribution, more states are occupied, 
and more energy eigenstates are visible in the spec- 
tral function (see Fig. [5J right panel). 
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Fig. 2: Spectral function for NN t in the broad-band limit D = 20t. Left panel: J = 0; the non- 

interacting levels are broadened either by the hybridization V or by a tiny cut-off 5/£ = 0.01. Middle 
panel: V = for different values of J and (3t = 20. One can clearly see that the gap A is controlled by t 
or J in the weak and strong coupling regimes, respectively. Right panel: Temperature dependence of the 
spectral function at J = 3t and V = 0. 



Moreover, due to the Hund's exchange that cou- 
ples locally the itinerant and localized spins degrees 
of freedom, non-local (short-range) magnetic cor- 
relations arise between the localized spins, so that 
they mutually influence their orientation. We will 
see that, depending on the ratio between J and 
V, and on the density of the itinerant electrons, 
those correlations can favor a FM or AF alignment 
of the localized spins. The exact summation tech- 
nique described in Sec. fallows us to calculate the 
probabilities P[{Si}] of Ea. (fT0| . in order to show 
which spin configurations are energetically more fa- 
vorable. As the probability of a configuration is -to 
a good approximation- only depending on the num- 
ber of FM and AF bonds, regardless of the exact 
relative position of these bonds, we can identify four 
representative spin configurations, for each of the 
hopping topologies introduced above. 

In the following we discuss general considera- 
tion in order to understand (or predict) what kind 
of magnetic correlation will affect the localized spins. 
In the isolated system (i.e. V/t = 0) only two en- 
ergy scales are competing, namely the hopping t 
which tends to minimize the energy and delocal- 
izes the itinerant electrons, and the Hund's cou- 
pling J, which couples their spins with the local- 
ized ones. Due to the Pauli principle, hopping pro- 
cesses between neighboring sites are only possible 
if the sites are occupied with one electron or two 
electrons with opposite spins. This leads to an ef- 
fective AF correlation between the itinerant elec- 
trons, and an AF alignment is particularly favored 
at half-filling, when there is on average one elec- 
tron per site. At the same time, in order to further 
minimize the energy, each localized spin tends to 
align to the spin of the itinerant electron (density) 
at the corresponding site. This mechanism leads to 
AF (superexchange) correlation between neighbor- 
ing localized spins. On the other hand, when the 
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Fig. 3: Probability P[{5j}] for four representative 
configurations (out of 2 6 ) of localized spins, classi- 
fied by the number of AF bonds, as a function of 
J and V, at fit = 20 and half-filling, in the NN t 
hopping topology. An AF kind of correlation of the 
localized spins is preferred for J > t and J > V. 



system is out of half-filling, the presence of empty 
(or doubly occupied) sites favors instead a FM dou- 
ble exchange correlation between localized spins. 
According to the considerations above, it is there- 
fore interesting to study in detail the NN t case 
at non- integer density, e.g., quarter-filling, as well. 
The role of the hybridization is to effectively de- 
couple the itinerant electron at each site from the 
other sites and the localized spin. The hybridiza- 
tion hence suppresses non-local magnetic correla- 
tions and leads to a more evenly distribution of 
spin configuration probabilities P[{Si}]. In the fol- 
lowing we analyze the probabilities of the different 
localized spins configurations as a function of V and 
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Fig. 

atures the system is characterized by more evenly 
distributed probabilities for the configurations of 
the localized spins. 
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Fig. 5: Same as Fig. [3j but at quarter-filling, where 
FM (short range) ordering is favored when J be- 
comes the dominating energy scale. 



J for NN t both at half- and quarter-filling. In Fig. 
|3]we show that for NN t at half-filling and f3t = 20, 
an AF (short range) ordering of the localized spins 
is energetically favored in a region where J > V, 
i.e. in the region where the probability of the con- 
figuration with six out of six bonds arranged AF is 
the highest (yellow region of the top left panel of 
Fig. [3]). When the temperature is raised to j3t = 5, 
states which are higher in energy become populated 
as well, leading to more evenly distributed proba- 
bilities, and therefore the AF correlation is weaker 
or confined to a smaller region of the V — J plane, 
as shown in Fig. 2J However, one also expects that, 
independently of temperature, the probabilities of 
the different configurations to become evenly dis- 
tributed even in the strong coupling limit / » (. 
This can be understood by considering that the 
short-range magnetic order arises due to the in- 
terplay of J and t. At J ^> t, the energy gain as- 
sociated to hopping processes can be estimated in 
second order perturbation theory as AE « — t 2 /2 J, 
and is vanishing for J/t —> oo. As a consequence, 
the localized spins become uncorrelatcd in this limit 
and non-local correlations are suppressed. 

In the quarter-filled case we find a similar be- 
havior as for the half-filling, but, as expected, with 
dominant FM rather than AF configuration, i.e. the 
highest probability corresponds to the configura- 
tion with zero AF bonds (bottom right panel of Fig. 
EJ. This result can be understood by considering 
that, at low densities, the itinerant electrons have 
an enhanced (with respect to half-filling) probabil- 
ity to hop to empty sites. These kind of processes 
are favored if the neighboring localized spins are 
aligned FM, and suppressed if they arc aligned AF, 
since the process would cost an interaction energy 
of order 2J. Recently the double exchange model 



for a finite system for fully saturated ti g spins has 
been analyzed in Ref. [14] . 

As we aim to compare the exact solution to 
DMFT, we also notice that DMFT has the prop- 
erty of becoming exact in the limit of infinite con- 
nectivity z — > oo [B]. One expects therefore that, 
increasing the connectivity of the system, non-local 
correlations to be averaged out and the mean field 
approximation to become reliable. As already men- 
tioned, a possibility to do this, without changing 
the geometry of the nanostructure, is to introduce 
longer range hopping. It is therefore interesting to 
study how the system changes when we go from 
NN t to all t, i.e., increasing the connectivity from 
z = 2 to z = 5. In Fig. [6] we show that indeed 
the probabilities of the different spin configurations 
are in general more evenly distributed than for the 
NN t case, even at low temperatures, while AF 
short range order is still slightly favored in the 
strong coupling limit. 



2.2 Comparison with DMFT 

Let us turn to the DMFT solution, which for the 
specific case of the double exchange model is tem- 
perature independent (see Eq. (fT4")) ) and generally 
neglects non-local correlations. We hence expect 
DMFT to describe a system better if there is no 
(short-range) ordering of the spins, i.e., when the 
different spin configurations are rather equally im- 
portant. 

For the NN t system, the left panel of Fig. [7] 
shows that the spectral function of the exact solu- 
tion exhibits many peaks at energies correspond- 
ing to the N2 N eigenvalues of the Hamiltonian, 
while DMFT treats these states within mean field, 
and therefore predicts broad bands. Except for this 
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Fig. 7: Representative results of the comparison of the spectral function A(u>) between the exact solution 
and DMFT. Left panel: J = 1.5t, V/t = 0, and /3t = 5 for NN t at half-filling which implies particle 
hole-symmetry for NN t. In the low- hybridization regime the exact solution displays many peaks which 
correspond to the eigenstates of the Hamiltonian, while DMFT rather predicts broad bands. Middle 
panel: J = 3i, V = t, and fit = 20 for NN t at half-filling; increasing the hybridization the DMFT 
description improves even at strong coupling. Right panel: J = t, V = t, and (3t = 20 for all t; increasing 
the connectivity the DMFT reproduces the exact solution, even in a non-perturbative regime. A similar 
behavior is observed also for quarter-filling (in the NN t topology) and in all other parameter regimes 
investigated. 
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Fig. 6: Same as Fig.[3]but in the all t hopping topol- 
ogy. Increasing the connectivity of the system leads 
to much more evenly distributed probabilities. 
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Fig. 8: Local (Su, left panels) and non-local (for 
neighboring sites Zi^+i, right panels) self-energy 
for J = 3i, V = t, fit = 20, and NN t hopping topol- 
ogy- 



striking difference, the spectral function is quali- 
tatively reproduced, and the DMFT spectrum is 
similar to a broadened version of the exact spec- 
trum. In the presence of the hybridization V, the 
peaks get broadened anyhow. Hence, the finite-1^ 
spectrum predicted by DMFT is more similar to 
the exact one, as can be observed in the middle 
panel of Fig. [7] In the all t topology, where non- 
local correlations are washed away due to the high 
connectivity, the exact spectral function is almost 
exactly reproduced by DMFT, as shown in the right 
panel of Fig. [7] 

A more precise characterization of the non-local 
spatial correlations between neighboring localized 
spins can be achieved considering the self-energies. 
In Fig. [5] and O we compare the exact and DMFT 
local (Sa) and the non-local (between neighbor- 



ing sites, Ei^i+i ) self-energies corresponding to the 
spectral functions shown in the middle and right 
panel of Fig. [71 respectively. 

In a Fermi liquid ground state, the local self- 
energy displays, around the Fermi energy Ep, a 
typical behavior Su(0) ~ —au> + 1(7 + Slu 2 ) with 
the coefficients a > and 7, i5 < 0. In the NN t 
case, shown in Fig. the DMFT local self-energy 
is highly non-Fermi liquid, displaying a large max- 
imum (in absolute value) in the imaginary part 
of £ and a positive slope of the real part (a < 
0). These are the fingerprints of the depletion at 
Ep observed in the corresponding spectral function 
A(lo), which is also a consequence of the band gap 
that characterizes the non-interacting ring. Inter- 
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Fig. 9: Local (Uu, left panels) and non-local (Uij, 
equivalent for each pair i^j, right panels) elements 
of the self-energies for J=t, V = t, and fit = 20, for 
both all t hopping topology. 



estingly, the exact local self-energy behaves, in con- 
trast, like in a Fermi liquid, and hence the compari- 
son with DMFT looks poor. In absence of non-local 
correlations, the exact local self-energy would cor- 
respond to a scenario in which A(lu) would display 
a (renormalizcd) quasi-particle peak. The exact so- 
lution also predicts a pseudogap. However, this is 
generated by non-local magnetic fluctuations, e.g., 
by Rei7i i j_)_i(0) between nearest neighbors (cf. also 
[T2irilll[r5] ). Non-local spatial correlations extend 
also to longer distances, and the corresponding el- 
ements of the self-energy (not shown) arc quan- 
titatively comparable to the nearest-neighbor one 
shown here. 

In the all t case, shown in Fig. [5J the picture is 
qualitatively the same, with the fundamental dif- 
ference, that the non-local contributions to the self- 
energy are sensibly smaller than the local ones, thus 
justifying that the DMFT spectral function well 
agrees with the exact one. Remarkably, due to the 
enhanced connectivity, the non-local self-energies 
remain small also at lower values of V/t, when the 
suppression of magnetic correlations due to the hy- 
bridization is less effective (not shown). This shows 
that for high enough connectivity, or strong hy- 
bridization, the Hund's coupling most likely does 
not result in a phase with short range order, and 
therefore can be treated reliably within a dynami- 
cal mean field approximation. 



3 Conclusion and outlook 

We have studied the double exchange model on a 
nanoscopic lattice. The advantage of this simple 
model is that it can be solved exactly also for finite 



systems. This allows us to systematically study and 
understand the nature of approximations such as 
DMFT. Hence we can address the question: Where 
is this approximation sufficient, where does it fail 
quantitatively or qualitatively? We found DMFT 
to be reliable in a wide range of parameters. In 
particular we investigated the role of non-local cor- 
relations which are neglected in any DMFT-likc 
calculation. These non-local magnetic correlations 
are suppressed if, e.g., the nano-cluster strongly hy- 
bridizes with metallic leads, or if it has a large coor- 
dination number, explaining the overall good qual- 
ity of the DMFT predictions. After our analysis we 
can conclude that DMFT can be a suitable tool to 
study also more complex correlated nanostructures 
in this parameter range. However, if non-local cor- 
relations are not negligible, DMFT is not reliable 
anymore, and a better approximation is needed. In 
this respect, it is possible to extend the method in 
the spirit of D/A [SUHUMIElinillll , which can ac- 
count for those non-local correlations that are due 
to local fully irreducible two-particle vertices. This 
step is expected to improve the approximation sig- 
nificantly. 
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